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Limitations of the modulation method to smooth a wire guide roughness 
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It was recently demonstrated that wire guide roughness can be suppressed by modulating the 
wire currents [H so that the atoms experience a time-averaged potential without roughness. We 
theoretically study the limitations of this technique. At low modulation frequency, we show that 
the longitudinal potential modulation produces a heating of the cloud and we compute the heating 
rate. We also give a quantum derivation of the rough conservative potential associated with the 
micro-motion of the atoms. At large modulation frequency, we compute the loss rate due to non 
adiabatic spin flip and show it presents resonnances at multiple modulation frequencies. These 
studies show that the modulation technique works for a wide range of experimental parameters. We 
also give conditions to realise radio-frequency evaporative cooling in such a modulated trap. 

PACS numbers: 39.35. +k,03.75.Be 
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I. INTRODUCTION 

Atom-chips are a very promising tool for cooling and 
manipulating cold atoms 01 • Diverse potentials, vary- 
ing on the micron-scale, can be realized and very high 
transverse confinements are possible. Envisioned appli- 
cations range from integrated guided atomic interferome- 
try 0,0; 01 to the study of low dimensional gases 
To take benefit of the chip technology, the atoms should 
be brought close to the current carrying wires. But the 
atoms then ex per ience a rough potential due to wire im- 
perfections [qI. I10I| and this used to constitute an impor- 
tant limitation of the atom-chip technology. However, a 
method to overcome this roughness problem, based on 
modulated currents, was recently demonstrated Due 
to the important envisioned applications of this method, 
a study of its limitations is crucial. 

The method to suppress atomic wire guide roughness 
relies on the proportionality between the longitudinal po- 
tential roughness and the wire current [ll|. Modulating 
the wire current around zero at a large modulation fre- 
quency, the atoms experience the time-averaged poten- 
tial, exempt from roughness. The modulation frequency 
uj must be large enough so that the atomic motion cannot 
follow the instantaneous potential. On the other hand, 
Lu should be small enough in order to prevent losses due 
to spin-fiip transitions. In this paper, we present analy- 
sis that go beyond the time-averaged potential approach 
and we identify the limitations of this method, both for 
low and large uj. We also investigate the possibility of 
using the radio-frequency evaporative cooling method in 
such a modulated trap. 

In sec. II, we present the considered situation. In 
sec. Ill, we investigate the Hmitations of the method that 
arise at small modulation frequency. We show that the 
atomic cloud is submitted both to a heating and to a 
residual rough potential and we quatitatively study these 
phenomena. In sec. IV, we compute the expected spin flip 
losses due to the time modulation of the magnetic field 
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FIG. 1: A wire guide produced by three current carrying 
wires. Geometrical deformation of the wire produce a lon- 
gitudinal potential roughness (of correlation length Ic) pro- 
portional to the wire current as depicted in the graph. 



orientation, that arise at large to. Finally, the last section 
gives some insights on the possibility of realizing radio- 
frequency evaporative cooling in the modulated guide. 



II. WIRE GUIDE 



A wire guide is obtained combining a transverse 
quadrupolar field and an homogeneous longitudinal mag- 
netic field Bq. The quadrupolar field can be reaHsed us- 
ing for example three current carrying wires as shown in 
FigdJ Because of wire deformations [lol. [l^. the current 
density inside the wires acquire non zero transverse com- 
ponents. This produces a longitudinal rough magnetic 
field bz proportional to the wire current, much smaller 
than the external field Bq. The method to effectively re- 
move roughness consists in modulating the wire currents 
at a frequency lu while the longitudinal field Bq is kept 
constant. 
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III. EFFECT OF THE MODULATION ON THE 
LONGITUDINAL MOTION 

Let us first assume the modulation frequency of the 
wire currents is small enough so that the atomic spin ori- 
entation can follow adiabatically the magnetic field ori- 
entation. The atoms are then subjected to the instanta- 
neous potential ^|B|, where /-i is the atomic dipole mo- 
ment. For Bq much larger than E/^ where E is the typ- 
ical transverse atomic energy, the instantaneous trans- 
verse potential is harmonic and proportional to the in- 
stantaneous wire currents. Since the oscillation frequency 
of this potential is modulated in time, the transverse clas- 
sical dynamics is described by a Mathieu equation, which 
predicts stable motion as long as w > 0.87u!±, where 
is the maximum instantaneous transverse oscillation fre- 
quency. This classical criteria is also predicted by quan- 
tum mechanics since the Wigner function evolves as a 
classical phase space distribution for a harmonic poten- 
tial In this paper, we assume this stability condition 
is fulfilled and we focus on the longitudinal motion, as- 
suming the longitudinal dynamics is decoupled from the 
transverse one. The longitudinal instantaneous potential 
is 

V{z,t) = u{z)cos{ujt), (1) 

where u{z) = fibz{z), sketched in FigdJ is produced by 
wire deformations. The idea of the method to smooth the 
roughness is that the longitudinal motion of the atoms 
does not have time to follow the time evolution of the 
potential. The atomic motion is then well described by 
the effect of the conservative potential {V{z,t)), where 
the time average is done over a modulation period. Since 
{V{z,t)) = 0, the atoms do not experience any rough- 
ness. We study below the conditions on the modulation 
frequency lu for such an approach to be valid. 

As the Hamiltonian experienced by the atoms is pe- 
riodic in time, we use the well known Floquet repre- 
sentation briefly presented below for the situation 
considered here. A new quantum number np is intro- 
duced, which gives a relative number of modulation en- 
ergy quanta. The Hamiltonian in this representation is 
time-independent and contains two contributions. The 
first one, 

oo 

i?o = (P^/i^m) + hLunF)\nF){nFl (2) 

riF — — OO 

does not couple different Floquet subspaces. The second 
one, 

oo 

Hi= u{z)/2{\nF){nF + l\ + \7iF + l){nF\), {3) 

np — — oo 

couples adjacent Floquet subspaces. If the state 
of the system in the Floquet representation is 
J2nF \'<Pnp){t)\nF), where \ilJnp)it) gives the state of the 



system in the manifold of Floquet number nF, then 
the state of the system in the bare representation is 
J2nF iV'iii?)!^^"'"^"*- Expectation values of observables 
contain cross terms involving different Floquet numbers. 
However, as long as evolution on time scales much larger 
than 1/uj is considered, such cross terms (interference 
terms) average to zero and the different Floquet states 
can be interpreted as physically different states. A given 
state has an infinite number of Floquet expansions. In 
particular, it is possible to assume that the initial state 
is in the Floquet manifold of Floquet number uf = 0. 

Let us consider a state \pa,nF = 0) of momentum 
Po in the Floquet manifold uf — 0. The modulated 
rough potential u is responsible for two different phe- 
nomena. First, it induces a change rate of the atomic 
energy. This irreversible evolution is due to the contin- 
uous nature of the rough potential Fourier decomposi- 
tion: the state \po,nF = 0) is coupled to a continuum 
of momentum states of the adjacent Floquet subspaces 
UF = ±1 and this coupling to a continuum induces a 
departure rate from the initial state, associated to a rate 
of kinetic energy change. Second, the modulated poten- 
tial is responsible for the well known adiabatic potential 
experienced by atoms in fast modulated fields [15|. We 
show that this adiabatic potential is due to processes of 
order two in u that couple the state \po,nF = 0) to the 
states \pi,nF — 0) via the virtually populated interme- 
diate states \q,nF = ±1). 

In the first sub-section, we investigate the first phe- 
nomenum and compute the associated heating rate for a 
cloud at thermal quasi-equilibrium. In the second sub- 
section, we derive the adiabatic potential experienced by 
the atoms. In both sections, we emphasize on the case 
where the potential roughness is that obtained at large 
distance from a flat wire having white noise border fluc- 
tuations. 



A. Heating of the atomic cloud 

Let us suppose the atom is initially in the state \pq, uf) 
of momentum po in the Floquet manifold uf ~ 0. As 
shown in FiglJl this state is coupled by u to the con- 
tinuum of momentum states in the Floquet manifold 
Uf ~ ±1, which leads to a decay of the initial state pop- 
ulation. The momentum of the final states that fulfill en- 
ergy conservation in the Floquet s ubspace tij? = — 1 are 
±?iq+ where 5+ = ■\/fcQ + 2mui/h, ko = po/h being the 
initial atomic wave vector. Decay towards these states in- 
volves the Fourier component ±q+ — fco of u and increases 
the kinetic energy of the atom by Tiui. If fcg > 2mu!/h, 
there exist states in the Floquet subspace uf — I that 
have the same energy as the initial state. Th e momentum 
of those final states are ±?ig_ where g_ = -y/fep — 2mw / h 
and decay towards these states decreases the kinetic en- 
ergy of the atom by hiu. A perturbative calculation, iden- 
tical to that used to derive Fermi Golden rule, gives an 
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FIG. 2: Transitions responsible for a heating of the atomic 
cloud. The parabolas give the energy -ffo, given in Eq. ([2]), 
versus the momentum p for different Floquet manifolds np. 
The state of momentum po in the Floquet manifold np = 
is coupled to different momentum states of the Floquet man- 
ifolds np = ±1 by the rough potential u{z). 



dE TTOjm 

dt ~ 2h' 

-e(|fco|- 



energy exchange rate 

-[( ^(-fco + q+) + S{-ko~q+))/q+ 
- y^2^^){S{-ko + q-) + S{-ko - q-))/q-] 

(4) 

where S{q) l/(27r) / e"?'^(u(0)u(2))dz is the spectral 
density of u, characterized by the correlation length Ic, 
and 6 (a;) is the Heaviside function that is zero for a; < 
and 1 for a; > 0. The derivation of Eq. ([4]) is detailed in 
the appendix A. As pointed out in the appendix, Eq. ([4]) 
is not valid for initial momentum very close to V2rnhuj. 
However, the range of kg for which the formula is not 
valid is in general very small and we ignore this in the 
following. 

Apart from the rms amplitude of the roughness, which 
accounts only for a multiplicative factor in the rate of 
the energy change, three energy parameters enter into 
account: E^j = Twj is the energy quantum corresponding 
to the modulation frequency, Em = ml^uj^ is about the 
kinetic energy of an atom that would travel over Ic during 
an oscillation period, and Ec = pl/{2m) is the atomic 
kinetic energy. In the following, we consider two different 
limits for which we give simplified expressions for the 
energy change rate: the classical limit and the quantum 
low energy limit. 

Let us first assume that E^/Ec ^ 1 and 

Ei^eU^ / Ec^'^ 1. We show below that this two con- 
ditions ensures the validity of the classical behavior. The 
first inequality ensures that one can expand the quantity 
g+ and q- in power of muj/ (fik^). The second inequality 
ensures that the spectral components S'(— fcg — q±) are 
negligible compared to the two other ones and that the 
latter are well approximated using a Taylor expansion of 
S. The energy change rate then writes 

- [2u;yvlS{u/vo) + S'iiu/vo)Luyvl] 7r/(2m), 

(5) 

where Up = hko/m is the atomic velocity. This energy 
change rate does not depends on h and is thus a classi- 



dt 



cal result. It is obtained through a classical calculation 
of kinetic energy change computed after expanding the 
atomic trajectory to second order in u. Note that, using 
the classical expression Ec = mvQ/2, the two conditions 
Eu,/Ea < 1 and E^eW^/eV^ < 1 are verified in the 
limit where fi goes to zero, as expected for the validity of 
classical physics. 

Let us now consider the limit E^/E^ ^ 1 and 

EcEm"^ / Et!"^ <C 1, that we denote the quantum low en- 
ergy limit. The first inequality ensures that the Heaviside 
function in Eq. (0]) is zero and that g+ can be replaced 
by V2?iw in the denominator. The second inequality en- 
sures that this replacement is also valid for the argument 
of the S function. Then the energy change rate given in 
Eq. ([4]) reduces to 

dE _ 



5(-fco + ^/2muj/n) + S'(fco + yj2mujlh) 

(6) 

This is a quantum result, sensitive to the fact that energy 
exchange between the atom and the oscillating potential 
involves the energy quanta huj. In the limit where E^ ^ 
El^/E^ {ko ^ l/'c), it converges towards a finite value 



dE 
'dt 



TT\^2mto 

~7W~ 



S{^j2muj/n), 



(7) 



that does not depends on the initial momentum hko/m. 

Let us now consider a cloud initially at thermal equi- 
librium with a velocity distribution n{vo). The heating 
rate, obtained after averaging Eq. ^ over n{vo), is 



dT ^ [°° , ,dE , 



(8) 



where ks is the Boltzmann factor. Although the heating 
rate depends on the precise shape of the spectral density 
S, some general properties can be derived. 

First, although the energy change rate may be negative 
for some velocities, we show below that dT/dt is always 
positive. For a longitudinally homogeneous gas, this pos- 
itivity ensures the increase of the entropy, as required by 
the second law of thermodynamics in the absence of heat 
exchange with the cloud and without gaining information 
on the system. To demonstrate that dT/dt > 0, we per- 
form a change of variables in the four integrals obtained 
by substituting Eq. ^ into Eq. ^ to find 

^-B^ = TTUj^/m/h 

C ^S{q) {n{oj/q - hq/{2m)) - n{u / q + hq/{2m))) 
+ f 5(g) {n{nq/{2m) - co/q) n{co/q + hq/{2m))) 

(9) 

where Qo = ^j2Miuj/fi. For a thermal equilibrium distri- 
bution, n(v) is a decreasing function of \v\. Furthermore, 
the spectral density is a positive function. We thus find 
that dT/dt > so that the effect of the potential rough- 
ness is always a heating of the cloud. Eq. ^ also shows 
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FIG. 3: Heating rate of a cloud as a function of its temper- 
ature, for a modulated rough potential whose spectral den- 
sity is given by Eq. I|lip . Graph (a): classical predictions 
(dashed lines) and asymptotic behavior at low temperature 
given by Eq. lfT2|l (solid line). Graph (b): exact result for 
[^/(27r) = 200h/{md^) (solid line) compared with the quan- 
tum low energy prediction (dotted line), the asymptotic pre- 
diction of Eq. I|13p (dashed dotted lines) and the classical 
result (dashes). 



that the heating rate goes to zero at very large temper- 
ature, since n is then about flat over the explored veloc- 
ities. 

Second, for large enough temperatures one expects to 
recover the classical result and the heating rate should 
not depend on h. Then the heating rate depends only on 
the four independent quantities (ii^), Em, oj and fcsT. 
Since (w^) enters only as a multiplicative factor in the 
heating rate, using dimensional analysis we show that 
kudT/dt is the product of {v?) / {mujl'^) and a function of 
kBT/Em- As a consequence, if the function giving the 
heating rate versus T is known for a given value of w and 
Emi then the heating rate is known for any value of T , 
Em. and oj. 

Finally, at low enough temperatures, the heating rate 
is well estimated by substituting Eq. ([6]) into Eq. ([8]). One 
expects that the heating rate converges towards Eq. ([7]) 
when temperature becomes much smaller than Tiw and 

In the following, we give quantitative results in the case 
of a potential roughness obtained at large distance d from 
a flat wire whose borders have white noise fluctuations 
of spectral density J/. In this condition, the spectral 
density of u is given by 



S{k) 



j^iM^k^K.ikdr 



47r2 



(10) 



where Ki is the modifled Bessel function of the flrst kind. 
The typical correlation length of u is the distance above 
the wire d so that Em = muo'^d^. The mean square of 
the rough potential is {u^) ~ 0.044(/i//o/)2 J//d^. In the 
following we use {v?) as a parameter instead of J/. The 
spectral density of u then writes 



where a ~ 23. 

We flrst study the heating rate predicted by classical 
physics. This classical heating rate is plotted in Figl3l[a) 
as a dashed line. The temperature and heating rate are 
scaled to Em and uj(u^)/Em respectively in the graph 
so that the curves corresponding to the classical predic- 
tions are independent of the problem parameters. We 
observe the expected decrease to zero of the heating rate 
at high temperatures. We also observe a rapid decrease 
of the heating rate as temperature decreases, for tem- 
peratures much smaller than Em- The maximum heat- 
ing rate is about 2.1{u'^) / {mivd'^) and is obtained for 
the temperature ksTM — 0.07Em- For d = 5 /xm and 
cij/(27r) = 50 kHz, which are parameters similar to that 
of the experiment presented in Tm = 1.8 mK. Typi- 
cal cold atoms temperatures are much smaller than this 
value and it is thus of experimental interest to investigate 
in more detail the regime T ^ Tm- 

The decrease of the heating rate for T ^ is ex- 
pected since then the atoms move on a distance much 
smaller than the correlation length of the rough poten- 
tial during a modulation period. The atoms are then 
locally subjected to an oscillating force almost indepen- 
dent of z and the atomic motion can be decomposed into 
a fast micro-motion in phase with the modulation and a 
slow motion. The micro-motion being almost in phase 
with the excitation force, barely no energy exchange be- 
tween the atom and the potential arises on a time scale 
larger than the modulation period. More quantitatively, 
we can derive an analytical expression of the heating rate 
in the regime where T ^ Tm, which shows the decrease 
of the heating rate as temperature decreases. For such 
low temperatures, as shown a posteriori below, wave vec- 
tor of S that contribute to the heating rate are much 
larger than l/d so that we can replace the Bessel function 
Ki(x) in Eq. ifTT|) by its asymptotic value at large x. We 
then flnd that the integrand in Eq. ([8]) is peaked around 
vo = i^^^iksT^EmY^^ /m and the Laplace method gives 
the following approximation of the heating rate: 



ksdT ^ ^ ^i^"^) ( '^'^ ^-3(i5„./(2/cBT))i/^ 



dt 



Er, 



\kBTj 



(12) 



S{k) ^ a{u^)d{kdfKi{kdf 



(11) 



where /? ~ 0.36. This asymptotic function is plotted in 
Figl3](a) (solid line). It coincides with the exact classical 
result within 20% as long as ksT < 0.002Em- The above 
expression of vq and Eq. ^ validate the expansion at 
large x of the Bessel function Ki{x) for T <C Tm- 

The limit of validity of the classical results described 
above is given by E^/Ec < 1 and E^eU^/eI'^ < 1, 
where Ec — touq, vq being the typical velocity involved 
in the heating process. Using the above value for vq, 
the condition of valididty of the classical regime reduces, 
for E^ < Em, to ksT » E^. For w/(27r) = 50 kHz, 
we find that the classical regime fails for temperatures 
T <ti 2 /iK. At lower temperatures, quantum analysis is 
required to estimate the heating rate. 

For the parameters chosen above (rf = 5 /im and 
uj/{2tt) = 50 kHz), Em/E^ = 10^* and the heating 
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rate is exponentially small at temperatures smaller than 
a micro-Kelvin when classical physics fails (the term 
g-3(£;,„/(2fcBT))i/3 ([12]) is 3 X 10-31 for T 1 ^K.). 

Thus, in order to investigate the heating rate beyond the 
classical approximation, we consider a different situation 
for which Em/E^ is only equal to 200. This would corre- 
spond, for the same distance d = 5 /im, to a modulation 
frequency of only 1 kHz. The exact heating rate, which is 
computed by substituting Eq. Q into Eq. ([HI, is plotted 
in Figl3](6). This calculation shows that the classical re- 
sult is valid up to a factor of 2 as long as ksT > 0.2E^. 
At smaller temperatures, the classical result underesti- 
mates the heating rate. At temperature much smaller 



than Eu. 



the heating rate is well approximated 



by the predictions in the low energy quantum limit where 
Eq. ^ is valid. This prediction is represented as a dot- 
ted line in the graph. At temperatures much smaller 
than El/Em (i.e. for fc^T < h^/{md'^)), the heating 
rate converges towards Eq. |[7]). Assuming E^/E^ ^ 1, 
then the expansion of S at large wave vector can be used 
and Eq. ([7]) gives 



kedT/dt = C(w^> /Hmujd^/nfe-'^V"''^'^''/'' (13) 

where C ~ 4.0. This asymptotic value is plotted in FiglS] 
as dashed-dotted lines. The heating rate is equal to this 
limit up to a factor of 2 as soon as ksT < 0.2E^/ Em- 

The heating of the atomic cloud can easily be made 
small enough experimentally to have no noticeable ef- 
fects. Let us for example consider the situation, sim- 
ilar to the experiment in [l|, where d = 5 fim and 
yjjv?) = 50 nK. If the modulation frequency is as low 
as 1 kHz, then the maximum heating rate is 3 ^K/s and 
is obtained for a temperature of 700 nK. Thus, for such 
a low modulation frequency, the heating may be a prob- 
lem in experiments. However, as soon as the modulation 
frequency is increased to 50 kHz, as in [l|, the maximum 
heating rate is only of 64 nK/s and is obtained at a large 
temperature of 1.8 mK. At lower temperature, the ex- 
ponential decrease of the heating rate shown in Eq. fT2|) 
rapidly decreases the heatig rate to completely negligible 
values. 



B. Effective remaining potential 

In this subsection, we show that Raman processes (of 
second order in u), in which adjacent Floquet states are 
virtually populated, are responsible for the effective po- 
tential 



Vad = {d.uY /{Amoj^). 



(14) 



This potential is a well known classical result [15| that 
corresponds to the kinetic energy of the micro-motion 
of the atom. At large oscillation frequency, the micro- 
motion has an amplitude ^ ~ — ^^^/(tow^) cos(a;i) much 
smaller than the correlation length of u. It is in phase 



np 




Q hko Tiki 



FIG. 4: Raman transitions responsible for the adiabatic po- 
tential of Eq. I|14p . The second order coupling between two 
momentum states of wave vector ko and fci is sketched. The 
coupling produced by two Fourier components of u of wave 
vector q and q = k\ — q — ko are represented as dashed and 
solid arrows respectively. 



with the excitation force and has a kinetic energy Vad- 
In this limit, since the micro-motion is in phase with the 
excitation force, the energy transfer between the atom 
and the potential, averaged over a modulation period, 
vanishes. Energy conservation then shows that the slow 
motion of the atom is subjected to an effective rough 
potential Vad- It is well known that this effective poten- 
tial, due to the fast atomic micro-motion, is responsible 
for the confinement in rapidly modulated Paul traps. A 
quantum derivation of Vad has already been done in [l3| 
using a secular approximation. Here we give an alterna- 
tive derivation based on the Floquet representation. 

Let us compute the effective coupling between the 
states j/co) and |fci) of momentum Tik^f and fiki respec- 
tively, both being in the Floquet subspace np = For 
this purpose, we first investigate the effect of a given pair 
of Fourier components of u that couple the two previous 
states. Their wave vectors are q and q' = ki — k^ — q. Four 
processes are involved in the effective coupHng between 
|fco) and \ki) , as sketched in Fig|4]and the effective cou- 
pling is the sum of the four amplitudes. The precise ef- 
fective coupling between two "ground" states coupled via 
an intermediate level has been investigated in The 
authors show that the effective coupling is 14^2 /A where 
Vi and V2 are the coupling to the intermediate state and 
A is the difference between the energy of the intermedi- 
ate state and the mean energy of the two "ground" states. 
Using this result, we find that the effective coupling as- 
sociated with each process is 



UgUg'/{A{n^(ko + Kf)/{2m) ±huj- Eo), (15) 



where k is q or q' depending on the process and Eq 



6 



h (fcp + kf)/{Am). Adding the four amplitudes, we find 



^eff 



h'^iko + q)^/m~2EQ 

4 \{h^{ko + q)^/{2m)~ Eo)^ -n^u}^ 
h\ko + q')ym~2Eo 

{h\ko + q'f/{2m) - Eaf - ^2^2 

(16) 

Assuming that the kinetic energies of the final, initial 
and intermediate states are all much smaller than Tllo, the 
denominator can be simplified to fi^uj'^ and we obtain 



2ma>2 



(17) 



Doing the sum over the pairs {q,q'), we find that the 
coupHng between the momentum states is that realised 
by the potential of Eq. (fT4|) . 

An alternative way to derive the adiabatic potential 
is to use a dressed Floquet representation. As shown 
in appendix B, calculations are more straightforward in 
this representation and no detailed knowledge on the ef- 
fective coupling corresponding to a transition through a 
virtually populated state is required. 

The residual roughness given in Eq. lfT4| constitutes 
a limitation of the modulation method. It scales as 
(u^^ /{ml'^iu'^) where Ic is the typical correlation length 
of u. Thus it is much smaller than the initial roughness 
amplitude as soon as Urms ^ mllu'^ . In the case where 
the roughness potential spectral density is that obtained 
at large distance d from a wire having white noise bor- 
der fiuctuations of spectral density J/ , we obtain a mean 
value 

(Feff) = 0M8j/j^^ = l.l{u')/imuj'd'). (18) 

If the wire edges deformations have a gaussian probability 
distribution, th en the roughne ss of the remaining poten- 
tial is simply ^(V;2j) - (y^ff)2 = ^/2{V^s). For d = 5 ^lm, 
uj/{2'k) = 50 kHz and {u'^)/kB = 500 nK, the roughness 
of the effective remaining potential is as small as 0.09 nK. 



IV. LOSSES DUE TO SPIN-FLIP TRANSITIONS 

All the previous analysis assume the atomic spin can 
adiabatically follow the direction of the instantaneous 
field when the current is modulated. In this section, we 
investigate the conditions on the modulation frequency 
for this adiabatic following condition to be valid. Non 
adiabaticity induces losses via spin-fiip transitions to the 
untrapped states. Intuitively, we expect that the losses 
are small for a modulation frequency much smalller than 
the Larmor frequency. In the following we compute this 
loss rate, following calculations made in [l3| for the DC- 
case. 

Let us consider a spin-one atom in the modulated guide 
described in section II with w » wj^. We choose a co- 
ordinate system whose origin is at the quadrupole field 



center and of axis x and y at 45° with the quadrupole 
axis as depicted in Fig[TJ We assume the atomic mag- 
netic moment is ^J/7i where J is the atomic spin angular 
momentum of components Jx,Jy and Jz along x,y and 
z respectively. We note | ± 1) and |0) the eigenstates of 
Jz of eigenvalues ±Ti and 0, respectively. 

The magnetic field direction depends on the position 
and on time. We apply a spatially and time dependent 
spin rotation TZ{t,x,y) so that 7?,|1) points along the in- 
stantaneous local magnetic field direction. In such a rep- 
resentation, the Hamiltonian is 



H = 'R 



1 PI 

2m 



dTZ^^ 

■n + ujz/n + ih — ; — 7^ 

dt 



where 



U = ^iBq + muj]_ cos^ {ijjt){x'^ + y^)/2. 



(19) 



(20) 



Here we assume mu!\{x^ + y"^) A'-Bo so that the har- 
monic approximation is valid. We also neglect the ef- 
fect gravity. This later assu mption is relevant as soon 
as g <C l±iLj'j^ where l± = y/fi/ {muj±) is the harmonic 
oscillator length. This condition ensures that the time- 
averaged potential is barely affected by the gravity and 
that the acceleration of the atoms over the spatial exten- 
sion of the trapped state has a neghgible effect. 

We choose the rotation 7?, as a product of a rotation 
along X and a rotation along y. To first order in x and 
y, n is 



11 = 1 + i9j:Jx/h + iOyJy/n. 



(21) 



where the angles of the rotations along x and y are 
9x = —b'xcos{ujt)/Bo and 6y = b' y cos{ujt) / Bq respec- 
tively. Here b' = uj±y/2mJ3o/Ji is the quadrupole gra- 
dient at maximum current. Calculation to first order in 
dx.Oy gives 



2 2 

2m 2m 



(22) 



where 



V, = ^^^^^^iv.-^Py)ml\ + h.c. (23) 

-DO 

and h.c. stands for hermitian conjugate. Here, we ignore 
the state | — 1), which is relevant for low enough coupHng 
(see below). The term 14, due to the fact that TZ de- 
pends on the position, is responsible for spin-fiip losses 
in time independent ioffe magnetic traps [l3|. Within 
the approximations made, the position dependence of TZ 
has no effect on the Hamiltonian within the spin state 
|1) manifold: the Coriolis coupling analysed in [20|, that 
corresponds to a rotation frequency proportional to b'"^ , 
is not seen in this calcuation. 
Similar calculations give 



^ dTZ ^ fiujb' . , , , , , 

ih — —Tl = — — sin(tj<)(x - z?/)|0)(l| 



dt 



\/2Bo 



(24) 
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This term, due to the time modulation of the local spin 
orientation, may also produce spin-flip losses in modu- 
lated traps. The condition uj uj± ensures that the 
term of Eq. (|24| has an effect much larger than the term 
of Eq. (|23| and we neglect the latter in the following. 

As in the previous section, we use the Floquet repre- 
sentation. The Hamiltonian H in the manifold of spin 
state |1) is decomposed into the term 

oo 

Ho= (//2TO+TOt^i(a;^+y^)/4-fnF?iw)|l,ni?)(l,7iF| 

(25) 

and the term 

H2 = mujl{x^+y^)/'il ^ |1, + 2) (1, -f- /i.e. J 

\nF— — oo / 

(26) 

that couples the Floquet np state to the Floquet states 
np ±2. Here, \l,np) is the state vector of an atom in 
the spin state |1) with the Floquet number np. The term 
H2 is due to the part of Eq. (flOl) that is proportional to 
coa{2u!t). Since we assumed uj ^ uj±, the effect of H2 is 
weak and can be treated perturbatively. 

We will compute the loss rate of an atom initially in 
the spin state |1) of Floquet number np = and in the 
ground state ^0 of Hq. The term of Eq. l(24|) couples 
this trapped state to the untrapped spin states |0) of 
Floquet numbers np = ±1. The energy spectrum of the 
spin state |0), which is unaffected by the magnetic field, 
is a continuum. Coupling to this continuum leads to a 
departure rate from the initial state, provided the Markov 
approximation is fullfiled This approximation also 

ensures that the state | — 1) can be neglected. We will 
show below that losses to the Floquet manifold np = +1 
are much larger than losses to the Floquet manifold np ^ 
— I. Thus, we consider in the following the final states in 
the manifold np = +1. Since TZ does not affect 

the longitudinal motion, we concentrate on the transverse 
degrees of freedom and normalise 0o to JJdxdy\(j)o\^ = 1. 

In addition, because ih^^^^TZ(j)o is, up to a phase factor 
e'**, invariant under rotation of angle 6 in the xy plane, 
the losses to spin states are isotropic in the xy plane. It 
is thus sufficient to compute the departure rate towards 
a plane wave traveUing in the x direction. The final state 
wave vector is 



2m(fiBQ + huj±/V2 — hiL!)/h, 



(27) 



and the Fermi Golden rule gives the departure rate 
,2 







V2Tr 



?ifc?e-^''''^//(™"-L). (28) 



The departure decreases exponentially with the bias field 
Bq, as for a usual time independent ioffe trap How- 
ever, in the modulated trap, an additional exponential 
factor in u!/uj± refiects the fact that the Floquet level is 



increased by one while the spin is fiipped. This transi- 
tion is associated with the "emission" of a quantum of 
energy Ticu, given to the oscillating magnetic field. Equa- 
tion l(28|) also shows that the departure rate goes to zero 
for a modulation frequency very close to the frequency 
^Bo/Ti + iLj±/\/2, i.e. for vanishing kf. This cancellation 
is due to the fact that the coupHng term of Eq. ((24l) is 
odd in X whereas the initial state is even and the final 
state, whose wave vector is vanishing, is fiat. The depar- 
ture rate towards the Floquet state np = —1 is identical 
to Eq. l(28l) . Lu being replaced by —lu. Since we assumed 
UJ ^ Lu±, the loss rate towards the Floquet state np = —1 
is negligible compared to Eq. l(28|) . 

The condition w 3> and Eq. l(28|) show that the 
loss rate is exponentially small when cu reaches uJi = 
{liBo/h + ujj_/^/2)/3, the value for which the initial state 
has the same energy as the untrapped state of Floquet 
number np = 3 and of vanishing momentum. For modu- 
lation frequencies smaller than uji , second order processes 
resonantly couple the initial state to the untrapped state 
|0) of Floquet number nj? = 3. In theses processes, repre- 
sented in FigEl the term H2 of Eq. l(26|) first transfers the 
atoms in the virtually populated intermediate trapped 
state |1) of Floquet number np = 2 before the term 
-TZ realizes the transfer to the untrapped state |0) 



dt 

of Floquet number np = 3. Although H2 is weak, the 
loss rate associated with the second order processes is 
much larger than the exponentially small Tq. 

More generally, for a given modulation frequency, 
losses are dominated by transitions towards untrapped 
states of Floquet number np — 2n + 1 where n — 
E{{fiBo/h + ujj_/y/2)/{2uj)-l/2), the function E{x) be- 
ing the integer part of x. Those transitions correspond to 
processes where the atoms is first brought to the interme- 
diate state |1) of Floquet number np = 2n hj n transi- 
tions produced by the term H2 and is then transferred to 
the untrapped state |0) of Floquet number np ~ 2n + 1 
by the term iMTZ'^/dm of Eq. Perturbation the- 

ory gives an effective coupHng between the state |1) of 
Floquet number np = and the states |0) of Floquet 
number np = 2n + 1 which is 



Un 



nuh'{x-iy) 2 ( 2 

-I = ymio I \x 

n\2V2Bo 



y')/{16hLu)r. (29) 



The eigenstates of Hq in the virtual intermediate states 
do not appear because, since we assumed ui ^ the 
energy difference between the intermediate states is neg- 
ligible and a resummation is possible. 

The departure rate from the initial state towards the 
Floquet state n_F = 2n -I- 1 is computed from Un using 
the Fermi Golden Rule. As for the calculation of Fq, it is 
sufficient to compute the departure rate in the x direction 
and we obtain 



(30) 



I// dxdy{x — iy)e^^f^4>Q{x,y){x'^ + y^)"| 



where k 



f 



2m{^Bo + hujj_/\/2- {2n+l)hu!)/h is 
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spin state |0) 




= 3 



FIG. 5: Transitions responsible for spin flip losses. For the 
spin state |1), the potential energy term of Eq. I|25p is repre- 
sented as well as the energy of the ground state in the np = 
manifold. For the spin state |0), we represented, for each 
Floquet manifold np, the whole energy spectrum, which is a 
semi-continuum starting at an energy npTiuj. The transitions 
induced by the term of Eq. I|24p are shown as dotted arrows 
whereas transitions due to the term H2 of Eq. l[26ll are shown 
as solid lines. The initial state is the spin state |1) of Floquet 
number np = 0. For the two final Floquet states np = 1 
and np = 3, only the dominant processes are sketched, whose 
amplitude are Uo and Ui respectively, where Un is given in 
Eq. l[29|l . In this picture, the odd Floquet state of |0) the 
closest to resonance corresponds to np = Z and losses are 
dominated by Fi, where r„ is given in Eq. ifSOjl . 



the wave vector of the final state. Using the gaussian ex- 
pression for the ground state 4>o(x, y) in Eq. l(30|) . we can 
show that r„ is the product of a polynomial in fc/ and 
of the exponential factor e.-^^k']/(niu>^) ^ r^^ie minima 
of the polynomial correspond to destructive interferences 
between the amplitude probabilities of paths having dif- 
ferent intermediate vibrational states. Since hkj/m is 
reduced by 4a; when increasing n by one and since we as- 
sumed uj ^ uj±, the exponential factor ensures, as stated 
above, that the total loss rate is dominated by the de- 
parture towards the highest Floquet subspace. 

Figure[6]gives the departure rate of the trapped ground 
state as a function of uj for fiBo = 50huj±. We observe 
several peaks that reflect the resonance behavior at in- 
teger values of (/i-Bo + '^±/V^ — w)/(2a;). The height of 
the resonances goes down with the integer n as expected 
since the order of the transition increases with n. We 
verify that the loss rate is dominated by the losses to- 
wards the Floquet state of highest odd Floquet number, 
as expected. Between two resonances, we observe the ex- 
pected exponential decrease of the loss rate. We observe 
a structure in the loss rate for losses to Floquet state 
larger than one, as expected. 

The lifetime of a thermal Maxwell-Boltzmann distri- 



o 
cq 

a. 




20 



30 



40 



50 



FIG. 6: Loss rate from the vibrational ground state versus the 
modulation frequency. The longitudinal magnetic field and 
the quadrupole gradient are chosen so that ^Bo/{tiuj±) = 50. 



bution is obtained after averaging the loss rate over the 
thermal distribution. In this calculation, (/jq in Eq. i[30|) 
is replaced by the eigenstate (j)i{x)(j)j{y), where i and j 
denotes the vibrational level. Neglecting changes of the 
final energy h'^k'j/{2m) across the thermal distribution, 
we can compute the departure rate Tq for a Maxwell- 
Boltzmann distribution. More precisely, writing Eq. l(30|l 
as a fourth integral and using properties of Wigner func- 
tions we find, for temperatures fcsT ^ TT'UJ±, 



^(pBa~nuj}/{kBT} 



(31) 



Because of the crude approximation that the final energy 
does not depend on the initial state, this result is only 
valid up to a factor of the order of unity. 

Experimentally, spin fiip losses can easily be avoided by 
properly choosing the modulation frequency. For exam- 
ple, let us assume the transverse oscillation frequency of 
the instantaneous trap at maximum current is u!±/{2'!t) = 
50 kHz and the longitudinal magnetic field fulfills Bq = 
50huj±/^. If /i = fiB where /is is the Bohr magneton, 
this corresponds to Bq = 1.8 G and the Larmor frequency 
fiBo/h is 2.5 MHz. In these conditions, the loss rate is 
dominated the term Tq of Eq. (|28)) as long as the os- 
cillation frequency is larger than 0.84 MHz and, in this 
frequency range, it is smaller than 1 s~^ as soon as w < 
2.2 MHz. For w < 0.84 MHz, losses become dominated 
by transitions towards states of higher Floquet number 
and the loss rate is peaked at modulation frequencies 
close to integer values of ^BQ/{2fiLo)+ujj[_/ {2^/2uj) — 1/2. 
In particular, loss rate goes up to about 25 s^^ for a 
modulation frequency close to 0.8 MHz. Thus, the vicin- 
ity of this resonance should be avoided experimentally. 
Resonances of higher order are less problematic since the 
maximum loss rate they induce is smaller than 0.1 s~^. 



V. RADIO-FREQUENCY EVAPORATION IN 
THE MODULATED GUIDE 

In this section we present general considerations on 
forced evaporation in a modulated guide. Since evapo- 
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rative cooling is most efficiently realised in a 3D trap, 
longitudinal confiment is required. A 3D trap can be 
obtained from the modulated guide of section II by ap- 
plying a z-dependent constant longitudinal field Bq{z). 
Here, we consider evaporation in the transverse plane 
(xy) at a given z position and denote as Bq the longi- 
tudinal magnetic field. For this purpose, in addition to 
the previous trapping potential, we apply a weak radio 
frequency magnetic field polarized in the x-direction, of 
frequency iujip and of amplitude Bjip . We consider here 
an atom of magnetic moment where J is the atomic 

spin angular mometum. 

Let us first give simple predictions, that rely only on 
the fact that, because of the modulation at 2uj of the 
trapping potential, the atomic Larmor frequency is mod- 
ulated in time. The modulation amplitude Sfl increases, 
in the transverse plane, with the distance r from the trap 
center according to S^l = pib''^r'^ /{AhBo). Considering 
only the internal atomic dynamics at a given position, 
the modulation of the Larmor frequency is equivalent, 
within the rotating wave approximation, to a frequency 
modulation of the radio-frequency field. In this picture, 
the radio frequency spectrum consists of a carrier at the 
frequency wrf and side-bands spaced by 2ijj, the relative 
amplitude of the n'^ sideband with respect to the carrier 
being Jn{fJ.{b' r)'^ / {SBohujj) , where J„ is the Bessel func- 
tion of the first kind. Thus, for a given frequency of the 
applied RF field, the coupling to the untrapped state is 
resonant for the positions r„ such that 



iURF = nBo/h + ^ib'^rl/{ABo) - 2nw, 



(32) 



where n is a integer. The coupling between the spin states 
close to a resonance position r„ is 



K = VoJn{^^{b'r)y{8BohoJ)) 



(33) 



where Vq is the coupling produced by the radio-frequency 
field in the absence of modulation. Such a shell structure 
of the spin fiip transition resonances is characteristic of 
AC magnetic traps. For example, the same behavior is 
expected in TOP traps [l^ , where the Larmor frequency 
is also modulated in time. 

In the following, we verify the statements made above 
by a more rigorous derivation. As in the previous section, 
we consider the representation in which the spin up state 
points along the local instantaneous magnetic field. The 
RF field produces a term in the Hamiltonian experienced 
by the atoms which is, to first order in the angles 0x = 
—b'xcos{uot)/BQ and 9y = b'y cos{ujt)/ Bq, 



Hrf = M-Brf cos(cjRFt) Jx 

-liB^-pb'v/Bo cos{uJnit)cos{ujTiFt)J^ 



(34) 



The right hand side is divided in two terms. The first 
term (first line) corresponds to the usual coupling be- 
tween the spin states in the presence of the RF field. 
The second term (second line) appears due to the time 
dependence of TZ. As in the previous section, in the fol- 
lowing we consider the case of a spin one state and we 



restrict ourself to the two spin states |1) and |0). We 
then have Hrfi = Hup^ + i/flF2i where 

Hrf, - MSRFCos(t^R,Ft)(|l)(0| + \0){1\)/V2 (35) 

and 

Hrf2 = - ^^^^ ^ cos{ujmt)cos{ujnFt)\l){l\. (36) 
Bo 

To analyze the effects of the RF field, we use the Flo- 
quet representation where two quantum number are used: 
the Floquet number np associated to the modulation fre- 
quency uj and A^RF , the number of radio-frequency pho- 
tons. We consider an atom in the trapped spin |1) state, 
with A^RF radio frequency (RF) photons. Because Brp is 
weak, we only consider transitions involving a single RF 
photon and we only consider transitions to the quasi- 
resonant states where the spin is and the number of 
RF photons is (A^rf — !)• 

Let us suppose the initial trapped state is in the 71^=0 
manifold. The term Hrf^ couples the initial state to the 
spin state in the manifold |A^rf — l,np = >. This 
transition is resonant for the position tq given by Eq. l(32|) 
and the coupHng to the spin state is ^Brf/^/2. The 
initial state can also be transferred to the spin state 
in the |A^rf — ^^np = ±2 > manifolds by higher order 
processes. These transitions, resonant for the position 
r±i given by Eq. (|32|) . can occur via two kind of pro- 
cesses, represented in Figl7| in the case where the final 
state lies in the manifold np = 2. In the first process, H2 
of Eq. (|26|) couples the initial state to the spin 1 state in 
the manifold \NBF,np = ±2), which is then transferred 
by the term Hrp-^ to the spin final state (process a). 
In the second kind of processes (process b), the trans- 
fer from the spin state 1 to the spin state is ensured 
by the term ih^^^^Tl of the Hamiltonian (see Eq. {Ml), 
and the term Hrp^ of the radio-frequency coupling is in- 
volved. In the case where to <^ ^Bo/h, the process (b) 
have a negligible amplitude and only the process (a) is 
important. 

In a more general way, the initial state can be trans- 
ferred to the final state of odd Floquet number np = 2n, 
the transitions being resonant at the positions r„ given 
by Eq. l(32|) . The dominant processes involve the first 
term Hrp given in Eq. (p2| and the term H2 to order n. 
The effective coupling between the trapped state and the 
spin state of Floquet number 2n, computed to lowest 
order in H2, is 



Vr 



n,eff 



(-1)" 



^ib'^r^ 
IGBnhuj 



\n\ 



HB 



RF 



V2 



(37) 



We recover here the result of Eq. I|33l) , in the limit con- 
sidered here where huj 3> ^6'^r^/(8i?o)- Thus the simple 
description in terms of frequency modulation of the Lar- 
mor frequency is sufficient to describe the physics. 

In conclusion, we have shown that the radio-frequency 
field is resonant for different trap locations, whose poten- 
tial energy differ by 2Tujj/kB. For a modulation frequency 
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Energy 



Processes (a) : 
l,nF = 2> 

l,nF= 1> 
Jl,nF = 0> 




Processes (b) 



|l,nF = 0> 



Nrf 



Nrf+1 



-|0,nF = 2> 



Nrf 




Nrf 



FIG. 7: Processes involved in the transition from the trapped 
state \l,nF — 0,Nrf} to the untrapped state \0,nF — 
2, Nrf — 1). The process (a) involves the Hamiltonian iif2 of 
Eq. JSSl) (solid line) and the term Hrf^ of Eq. (dashed 
line). The process (b) involves the term Hrf^ of Eq. II36P 
(dashed line) and the term ihdK-^/dtTZ of Eq. ^ (solid 
line). For uj ^ /iBo, the process (a) is the dominant one. 



of 50 kHz, the potential energy difference between two 
resonances is 3 fiK. For a temperature of the order of 
3 fiK or higher, some resonances are present inside the 
atomic cloud and induce spin-flip losses. To overcome 
this problem, a precooling stage in a static trap down to 
temperatures smaller than 3 /iK is required. For clouds 
whose temperature is smaller than 3 /iK, the evaporation 
process in the modulated guide involves only one radio- 
frequency knife, so that the evaporative cooling is similar 
to that realised in a trap made by DC currents. Choos- 
ing the frequency of the radio frequency field so that the 
transition involved in the cooling process is the transition 
that does not change the Floquet number is interesting 
for two reasons. First, as shown in Eq. (|37l) . the cou- 
pling between the trapped and the untrapped state of 
this transition is larger than that of higher order tran- 
sitions that change the Floquet number. Second, this 
coupHng is homogeneous and is thus constant when ujrf 
is chirped. 



VI. CONCLUSION 

The careful study of the limitations of the modulation 
technique to smooth wire guide roughness performed in 
this article show that this technique is very robust and 
accepts a wide range of modulation frequencies. More 
precisely, on one side, we have shown that the unwanted 
effects of the modulation on the longitudinal motion are 
neghgible for realistic parameters as soon as modulation 
frequency is larger than 10 kHz: both the heating of 
the cloud and the remaining effective roughness are very 
small. On the other side, the calculation of the losses due 
to spin flip transitions shows that, for reahstic parame- 
ters, these losses are neghgible as soon as the modulation 
frequency is smaller than a few hundred of kHz. 

The modulation technique is thus a very promis- 
ing tools that should enable to fully take advantage of 
atom chips devices. In particular, the study of one- 
dimensionnal gases in the strong interaction regime on 



an atom chip can be considered. 
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APPENDIX A: DERIVATION OF THE ENERGY 
CHANGE RATE 

The derivation of the heating rate follows that of the 
Fermi Golden rule. For the calculation, we assume a 
quantification box of size L and periodic boundary con- 
ditions. We assume the atom is initially in the state of 
momentum po in the Floquet subspace np = 0. For sim- 
plicity we consider only the transitions towards the Flo- 
quet state np = —1. After a time t much smaller than 
the departure rate, the change in kinetic energy AE can 
be deduced from perturbation theory and we obtain 



(Al) 



Here Uq = J dzu{z)e'^'^^ / L is the Fourier component of 
wave vector q of u{z) and 



f[q, i) = (e - 2 ^ 



(A2) 



where e = huj + h^q^ /{2m) + hpQq/m is the energy change 
associated to the transition involving the Fourier com- 
ponent of u of wave vector q. The terms Uq are com- 
plex random numbers without correlation between them 
and of mean square value = S{q)L/{2TT) where 

S{q) = l/(27r) / dze''^^u{z)u{0)) is the spectral density 
of u. For a large enough quantification box, the term 
f{q) barely changes between adjacent Fourier compo- 
nents and one can replace J2q Wql^fi^) by / d.qS{q)f{q) 
in Eq. IjAip . For a time t large enough so that the 
function (e — huj)S{q{e)) is about constant on the in- 
terval e e [—h/t,h/t], the term sin(et/2)/e^ can be re- 
placed by the distribution t(5(e)7r/2. We then recover 
the first term of Eq. |[4]). The previous condition on t 
and the condition that t is much smaller than the de- 
parture rate F can be fulfilled simultaneously only if the 
function (e — huj)S{q{e)) is about constant on the inter- 
val e e [— r,r]. This is the condition of the Markovian 
approximation. This condition is fulfilled provided that 
both S{q) and its correlation length are small enough. 

The calculation is similar for losses towards the Flo- 
quet manifold np — 1 and one finally recover Eq. 
The Markovian condition for the transition towards the 
Floquet manifold np — 1 is n ot fulfilled for initial mo- 
mentum very close to V^mfiuj since the atoms are then 
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sensitive to the fact that the continuum is not infinite for 
e < 0. Similar non Markovian situations have been stud- 
ied, for example, in photonic band gap materials j2l| and 
oscillatory behavior and decay towards a non vanishing 
population of the initial state are expected. 

APPENDIX B: ADIABATIC POTENTIAL IN 
THE DRESSED STATE REPRESENATION 

In this appendix, we rederive the adiabatic potential 
given Eq. (fT4| using a dressed representation, where a 
local z-dependent unitary transformation 0{z) is applied 
to the Floquet states so that the resulting dressed states 
\n){z) are eigenstates of the potential energy part of the 
Floquet Hamiltonian (term hujn p of Hq given in Eq. ([2]) 
and term H2 of Eq. ([3])). By symmetry, the energy of 
the dressed states \n){z) is nfiuj^ as that of the bare Flo- 
quet states. Using the properties of the Bessel functions 
{Jk+i{x) + Jk-i{x))x/{2k) = Jkix) and ^„ Jn{x)^ = 1, 
we show that the decomposition of \'n){z) in the un- 
dressed Floquet basis {\k)o) is 

00 

|n)(z)= Jk{u{z)/u;)\n + k)o^O{z)\n)„. (Bl) 

k——oc 

This well known result has been used in several other 
situations jl^l- In the dressed state representation, the 
state of the system is ip ~ O^^ipo where ipo is the state 
of the system in the undressed representation and the 



momentum operator, p = O ^pO, writes 

p^p-Y,{k,z\ind,\n,z)\k){n\ (B2) 

n,k 

where p ~ —ikd^ is the momentum operator that pre- 
serves the Floquet number and dz is a short notation 
for d/dz. Thus, in the dressed state representation, the 
Hamiltonian is decomposed into three terms: 

Ho ^ p'^ /2m ~ nphuj (B3) 

that does not couple different Floquet states, 

H, = -h/{2m) (pE„„„. \ni){n,\zdz\n2){n2\ (B4) 
+ E„i,„2 \ni){ni\idz\n2){n2\p'j , 

and 

H2=-h?/{2m) \ni){n2\{n^\dz\n^){n^\dz\n2) . 

ni ,712 ,713 

(B5) 

Since = {Jk-i - Jfe+i)/2 and Y^k^^Jk+n = ^n, Hi 
couples adjacent Floquet states. This term is responsi- 
ble for the heating of the cloud. On the other hand, H2 
contain a term H2,ad that does not change the Floquet 
number. Using the above properties of the Bessel func- 
tion, we find that H2.ad is just the adiabatic potential of 
Eq. in]). 
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